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The tunneling current between independently contacted graphene sheets separated by boron ni- 
tride insulator is calculated. Both dissipative tunneling transitions, with momentum transfer due 
to disorder scattering, and non-dissipative regime of tunneling, which appears due to intersection 
of electron and hole branches of energy spectrum, are described. Dependencies of tunneling current 
on concentrations in top and bottom graphene layers, which are governed by the voltages applied 
through independent contacts and gates, are considered for the back- and double-gated structures. 
The current-voltage characteristics of the back-gated structure are in agreement with the recent 
experiment [Science 335, 947 (2012)]. For the double-gated structures, the resonant dissipative tun- 
neling causes a ten times enhancement of response which is important for transistor applications. 
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I. INTRODUCTION 

In contrast to the tunneling processes between bulk 
materials, [Ij] the tunneling between low-dimensional sys- 
tems must be assisted by scattering in order to satisfy the 
momentum and energy conservation laws, see results and 
discussions for double quantum wells or wires in Refs. 2 
or 3, respectively. When the splitting energy between 2D 
states (this energy A is determined by transverse voltages 
applied across structure) exceeds the collision broadening 
energy F {Ti/T is the departure time), the tunneling prob- 
ability appears to be proportional F. In conditions of tun- 
neling resonanse^ when A <^ F,this probability is propor- 
tional to fi/T, [4] i.e. the tunnel current depends on the 
scattering time in the same way as the current in metallic 
conductor. The breakdown of the dissipative tunneling 
regime is possible if the energy spectrum branches are 
intersected and the energy-momentum conservation laws 
are satisfied without scattering. For example, the in- 
tersection of the parabolic electron branches in double 
quantum wells takes place if the magnetic field is applied 
perpendicular to the tunneling direction, see Q and Q 
for the experimental data and theory. Similar intersec- 
tion between the linear branches of gapless energy spectra 
should take place in graphene/boron nitride/graphene 
(G/BN/G) heterostructure. Such a structure was re- 
ported recently [7|] and the tunneling transistor, which 
is based on the independently-contacted G-sheets con- 
nected through a few monolayer BN, was demonstrated. 
Q In contrast to the semiconductor heterostructure case, 
the independently-contacted G/BN/G structures can be 
easily realized with the use of the single-layer-transfer 
technology. But a complete theoretical investigation 
of tunneling current in such a structure is not performed 
yet (some numerical results on the tunneling conductance 
were reported recently but the current-voltage char- 
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FIG. 1: (Color online) (a) Gt/BN/Gb structure under volt- 
ages Vtb and Vhg applied to top (t-) and bottom (6-) layers 
through independent contacts; back gate (blue) is separated 
by substrate of thickness ds and top gate under voltage Vtg 
is shown by blue dashed lines, (b) Dispersion laws of t- and 
b-layers (black and gray crossed lines) with Fermi energies spt 
and epb for the electron-electron tunneling regime, (c) The 
same as in (b) for the electron-hole tunneling regime. 

acteristics were not analyzed) and a problem is timely 
now. 

In the paper, we calculate the tunneling current / be- 
tween the independently-contacted top (Gt) and bottom 
(Gfc) graphene layers separated by BN. We analyze the 
dependencies of / on the sheet concentrations (Fermi en- 
ergies) and on A, which are determined by the gate volt- 
ages applied to the contacts, see Fig. la. Depending on 
voltages applied, one can realized either electron-electron 
(hole-hole) or electron-hole regimes of tunneling, as it is 
shown in Figs, lb and Ic, respectively. In the latter 
case the cross-point i? = is located between the Fermi 
energies {ept > > epb or vise versa) and the non- 
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dissipative regime of tunneling takes place in addition 
to the resonant dissipative tunneling transitions. The 
current- voltage characteristics appears to be different for 
these regimes. For the back-gated structure, the results 
are in agreement with the experimental data. Q For 
double-gated structures, the resonant dissipative tunnel- 
ing regime can be realized, with a ten times enhancement 
of response. 

The paper is organized as follows. In Sec. II, 
we present the basic equations which describe the two 
regimes of interlayer tunneling. In Sec. Ill we analyze 
the current-voltage characteristics and compare the re- 
sults for the back-gated structure with the of experimen- 
tal data. [1] The last section includes the discussion ap- 
proximations used, and the conclusions. In Appendix, we 
evaluate the effective tunneling Hamiltonian for G/BN/G 
structure. 



Similarly to Ref. 6, one express {p} through the carrier 
distributions determined by [p] and the tunneling current 
density takes form 
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dtSp ([p] 
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Here Sp(. . .) means both summations over states of car- 
riers in Gt- and G{,-layers and averaging over lateral 
disorder which should be included in the Hamiltonians 
hGt,b- Calculations of Sp(. . .) are performed below with 
the use of the basis formed by 2-row wave functions in 
fcth layer 4'x determined by the eigenvalue problem 
^Gfc^x'"'' — Sa'^x"'^- We introduce the spectral density 
matrix, labeled by I, I' =1 and 2, as 



II. BASIC EQUATIONS 
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Under consideration of Gt/BN/G;, structure, we use 
the tunneling Hamiltonian which connects Gt and Gf, 
layers described by 2x2 matrices hctfi, i-e. we introduce 
4x4 matrix [ill 
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Here we have separated the Hamiltonian of uncoupled 
layers, Hg, and the tunneling contribution, T, written 
through the 2x2 matrix f = f + which is determined by 
a stacking geometry of the structure, see Appendix. The 
charge density in Gt and Gf, layers, q[''^ (here and below 
k = t,b), and the tunnel current density. It , are deter- 
mine through the 4x4 density matrix pt by the formulas 
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Here is the normalization area, Pk is the projection 
operator on the fc-states, and the interlayer current op- 
erator J is determined from the charge conservation re- 
quirement It = dQf^ jdt = —dQ^^/dt (see similar calcu- 
lations in Refs. 6), so that 
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As a result, tunneling processes are described by the 
above-introduced matrices T and /, as well as the density 
matrix governed by the standard equation: ih{dpt/dt) — 

[Hgbng, Pt]. M 

Further, we separate the diagonal and non-diagonal 
parts of the density matrix pt = [pt] + {pt} which describe 
the distribution of carriers in Gt- and Gf,-layers and the 

tunneling current. It = (4/ L'^)Sp (^I{pt}^ ^ respectively. 



where the averaging over random disorder in fcth layer 
(. . .)fe is performed. Using the Fermi distribution for 
heavily-doped layers, when [p] is replaced by the 9- 
function 0{£Fk — £ka) with the Fermi energies Spk in 
Eq.(4), we transform / into 
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xtr f+ig'(x'-x)fi^^(x-x') 



where tr(. . .) means summation over the matrix variable. 
Below, we express in the momentum representa- 

tion as = « {gZ' - ^) /2-, where 

is the retarded Green's function of kth layer with the 
cross point energies ±A/2 written through the level- 
splitting energy. Within the model of short-range dis- 
order with the same statistically independent character- 
istics for k ~ t and b, the Green's function in the Born 
approximation takes form [l3| 



r<Rik) _ p(+)^(fe) I p(-)r^(fe) = n 



E.p 
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Ep 
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Ggp « [(E T A/2) (1 + As^A/2 + ^g) -vp] \ (7) 

where —A/2 and -l-A/2 are correspondent to t- and b- 
layers and v ~ 10® cm/s is the carrier velocity. The 
projection operators on the conduction (-I-) and valence 



(— ) band states, Pj 



(±) 



[1 ± (o- • p)/p] /2, are writ- 



ten through the 2x2 isospin Pauli matrix &. The self- 
energy contributions {E =F A/2)(A£;:p^/2 + *5) are writ- 
ten through the logarithmically-divergent real correction 
which is proportional to A^ = (f;7r) In (£'c/|i?|) and the 
coupling constant g. This approach corresponds to the 
short-range scattering with the cut-off energy Ec- 
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The interlayer tunneling is described by the parameters 
= tr(f+f), r2 = 1 ^tr(f+a^f<7j, (8) 



which appear under calculation of the matrix trace in Eq. 
(6). Using Eqs. (7) and (8) we transform the tunneling 
current density (6) into the form 



p 

+T,2lm^^pImeBp 
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(9) 



where the explicit expressions for Q^^^ and Q^^^ are de- 
termined by Eq. (7). Due to the in-plane symmetry of 
the problem, / is written as the double-integral over the 
p-plane and the energy interval {ept^^Fb)- 

Integrations in (9) are performed analytically for the 
coUisionless case, g ^ 0, and the tunneling current den- 
sity is given by the sum of the dissipative current, which is 
oc i5r(A), and the non-dissipative contribution (cx xl^D- 
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Here the resonant dissipative contribution (at ept^Fb > 
0) is written through the (5- like function (5r(A) — 
r/[7r(A2 + r^)] with the phenomenological broadening 
r ~ g{£Ft + £Fb)- At EFtSFb < 0, the non-dissipative 
contribution is written through the factor % (0 < x < 1 
because > T^) determined by the tunneling energies 
T and Tg. These parameters depend on the stacking 
order of the Gt/BN/Gf, structure and rough estimate 
for the case of A^-layer BN barrier [15|| is determined 
by Eq. (A. 4), see Appendix. As a result, we obtain 
T,Ts ^ 7(7/esjv)^, where 7 ~0.4 eV is the interlayer 
overlap integral and Ebn <2.5 eV is of the order of the 
c- and w-band energies in BN. We use below T,Ts ^13 
/ieV for iV = 4 and ~6 ^eV for 6-layer BN barrier, 
which are in agreement with the experimental data of 
Ref. 8. The scattering parameters used in Eqs. (9, 10) 
are taken from the conductivity measurements, see sim- 
ilar considerations in Refs. 16. We choose the cut-off 
energy Ec ^0.2 eV and the coupling constants g ~0.3 or 
0.15 corresponded to the maximal sheet resistance ~4 or 
2 kJ7 per square. 

The current density / is dependent on spt, £Fb and A, 
which are determined by the drops of voltages applied to 
three- or four-terminal structures, see Fig. la. Before 
study of the current-voltage characteristics, we consider 
the dependencies of the total current / on these three 
energies. In Fig. 2 we plot the tunneling current, which 
is determined by Eqs. (9) and (7), versus A for different 
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FIG. 2: (Color online) (a) Tunneling current density / versus 
level splitting A for the resonant tunneling regime at £f =200 
meV and Se =50 meV (1) or 100 meV (2) for g = 0.3. Solid 
and dashed curves are plotted for x =1 ^ind x =0, respec- 
tively. Curve (3) is plotted for parameters (1) at g — 0.15. 
Dotted curve fits dependencies (2) with the use Eq. (10). (b) 
The same as in panel (a) for the non-dissipative regime, if 
ef =0 and Se =200 meV (1) or 100 meV (3). Curve (2) cor- 
responds ef =200 meV, 5e =100 meV and dotted line gives 
oc A contribution in Eq. (10). 



Fermi energies determined by £f = (^Ft + £f&)/2 and 
Ssf — £Ft — £Fb at X ~1 01' 0- The resonant dissipa- 
tive regime of tunneling is realized at \eF\ > \deF/2\ and 
peak value /(A = 0) increases both with doping levels 
and with g, as it is shown in Fig. 2a. If \ep\ < \Sep/2\, 
the electron-hole tunneling regime takes place and the 
non-dissipative contribution becomes dominant with in- 
creasing of A, where / cx x|A|. At small A the de- 
pendency /(A) is transformed into narrow peak due to 
the dissipative contribution with weak broadening, see 
Eq. (10) and Fig. 2b. The numerical results given by 
Eq. (9) are in a good agreement with the approximations 
(10), shown by the dotted asymptotics, because of weak 
(<10%) contributions from the renormalization of energy 
spectra. 

Neglecting the quantum capacitance contributions and 
the near-contact drops of potentials, we use ept — £Fb — 
eVtb and the Gauss theorem connected the carrier concen- 
trations in graphene with interlayer electric fields, Vtb/d 
and Vbg/dg (here d and ds are thicknesses of BN layer 
and Si02 substrate). As a result, the Fermi energies and 
A are connected with drops of voltages as follows: 

2£Ft,6 ±eVtb - F{eVtbed) - F{eVbged, - eVtbSd), 
A ~ eVtb + F{eVtbed) - F{eVbged, - eVtbEd), 
— e{nvf 



F{x) — sign{x)^J\x\ 



Ae^d 



(11) 



Here and "— " stand for Gf- and G^-layers, eptfi > 
or < correspond to electron or hole doping, and the 
dielectric permittivity e ~4 is the same for BN and Si02 
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FIG. 3: (Color online) Tunneling current I versus Vtb at 
Vbg =Q for structures with 4- and 6-layer BN barriers (black 
and red curves, respectively; squares and triangles are exper- 
imental points [§]). Solid, dotted, and dashed curves corre- 
spond to X =1. 0.5, and 0. Insets show A and concentrations 
nt or rib (solid or dashed lines) versus Vtb- 




layers, see Refs. 7 and 8. For the double-gated struc- 
ture with top gate separated by BN insulator of thickness 
dt, one obtains similar expressions for ept.b and A after 
the replacement eVtbSd eVtbSd — eVtgSdt , see voltages 
shown in Fig. la. The approach (11) is valid for the 
heavily-doped graphene layers, so that the fields Vtb/d, 
Vbg/ds, and Vtg/dt should be strong enough. On the 
other hand, these fields are restricted by the breakdown 
condition for BN layer, when these fields should be ^ 7 
MV/cm. [13 



III. CURRENT- VOLTAGE CHARACTERISTICS 

In this section we analyze the current- voltage charac- 
teristics for the back- and double-gated structures. Be- 
low, the tunneling current density is determined by Eqs. 
(7) and (9) with the Fermi energies and the level split- 
ting written through the voltages applied according to 
Eqs. (11). The parameters described both the elastic 
scattering processes and the interlayer tunnel coupling 
are chosen the same as for the above calculations shown 
in Fig. 2. 

For the back-gated structure with A^-layer BN barriers 
[N —4 and 6) and g =0.3 the dependencies of / versus 
Vfb at fixed Vtg are shown in Fig. 3 at Vtg = 0, when 
I-Vtb — ~^Vtb- [3 The electron (hole) concentrations 
are oc Vtb and A is the sum of linear and square-root 
functions, see insets. The current-voltage characteristics 
liVtb) are in reasonable agreement with the experimen- 
tal data of Ref. 8 if we used the scale factor Jt with the 
above-estimated T and Tg. The dependencies on x are 
weak enough (<25% for N ==4 and <10% for N =6, not 
shown). If A ^0, when Vtb — -0.55 V for 6-layer barrier, 
a negative resonant contribution due to dissipative tun- 
neling gives '^30% variations of the I-V characteristic. 
Such a peculiarity was not found in Ref. 8; probably, it is 
due to a lateral redistributions of charges or a large-scale 
inhomogeneities of the samples used. The I-V charac- 
teristics of back-gated structures at Vbg ^ are plotted 



FIG. 4: (Color online) (a) Current density / versus Vtb for 
4-layer BN barrier at Vbg =0 (red), -25 V (blue), and -50 V 
(black) . Experimental points are shown for Vbg =0 (triangles) 
and -50 V (hexagons). Inset shows A versus Vtb. (b) The same 
as in panel (a) for 6-layer BN barrier. 



in Fig. 4 for N =4, 6 and x = 1, 0. Once again, there 
is a reasonable agreement of /(Vth,Vhg) with the avail- 
able experimental data for iV =4 at Vtb < 0. A visible 
asymmetry of / takes place if Vtb ~^ —Vtb but deviations 
from experimental data increase. The dependencies for 
N = 6 are similar but / is '-^2-times weaker and the res- 
onant dissipative contribution is absent for Vbg — —25 V 
because A 7^ 0, see inset. 

An enhancement of the resonant dissipative tunnel- 
ing contributions takes place in the double-gated (four- 
terminal) structure, when / depends on Vtb, Vtg, and Vbg- 
These dependencies are shown in Fig. 5 for the struc- 
ture with 10-layer BN cover layer (dt =3.4 nm; with a 
negligible tunneling, if dt ^2 nm [l^) at different top- 
and back-gate voltages {Vtg = corresponds to the three- 
terminal structure) . Since the condition A = is realized 
now at higher energies, the resonant dissipative tunnel- 
ing peaks are >10 times greater than the background cur- 
rent. This is the central result which should be important 
for the transistor applications. Note, that for A ~0 at 
Vtb "^0 (the curve 3 in Fig. 5a) the resonant dissipative 
peak is suppressed. Beyond of this narrow region, the 
resonant condition A(Vth) ~ should not be suppressed 
by lateral inhomogeneities and the peaks caused by the 
resonant dissipative tunneling should not be smeared. 

Overall, the consideration presented here gives an ad- 
equate theoretical description of the tunneling transistor 
effect which is in reasonable agreement with the exper- 
imental data. ^] The analysis performed opens a way 
for the verification of scattering mechanisms and tunnel- 
ing parameters in G/BN/G structures. In addition, the 
double-gated (four-terminal) structure was not consid- 
ered before and this structure show a great (ten times) 
enhancement of tunneling current tunability. 
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FIG. 5: (Color online) (a) Dependencies liVtb) for double- 
gated structure with 4-layer BN barrier and 10-layer BN cover 
layer under Vbg = -50 V and Vtg = 1.5 V (1), 0.75 V (2), - 
0.75 V (3), and -1.5 V (4). Inset shows A versus Vtb- Solid 
and dashed curves correspond to x =0 ^-nd li respectively, 
(b) The same as in upper panel for Vbg = V. 



IV. CONCLUSION 

We have adopted the theory of low-dimensional tun- 
neling 0, II] to the case of the tunneling charge trans- 
fer in independently contacted graphene structures with 
multi-layer BN barrier taking into account the two ba- 
sic mechanisms of charge transfer: resonant dissipative 
tunneling and non-dissipative tunneling. Both the non- 
dissipative tunneling current and the non-resonant dissi- 
pative tunneling processes are responsible for the current- 
voltage characteristics of back-gated structures. The 
resonant dissipative tunneling regime is achieved for the 
double-gated structures and these devices demonstrate 
a ten times enhancement of I-V characteristics which is 
important for transistor applications. 

Further, we discuss the assumptions used. Because of 
a luck of data on stacking order in G/BN/G structures 
the tunneling energies T and Tg were estimated from the 
current- voltage characteristics of Ref. 8 and this result 
is in agreement with the tight-binding model, [l^ More 
accurate estimates for T and should be based on ad- 
ditional structural measurements. The simplified elec- 
trostatics description, given by Eq. (11), fails for the 
low-doped G(- or G^-layers, under weak interlayer fields 
applied. In these narrow regions, a more complicated de- 
scription, which involves the quantum capacitance effect 
and the contact phenomena, should be applied. The rest 
of assumptions (the model of elastic scattering, [l^ iden- 
tical scattering in Gt^b-layers, weakness of long-range dis- 
order, and the single-particle approach) are rather stan- 
dard. 

To conclude, we believe that the description of tunnel- 
ing processes is an essential part of physics of graphene 



and the results obtained can be applied for characteriza- 
tion of scattering mechanisms and tunneling parameters 
in the tunnel-coupled graphene structures. More impor- 
tantly, that these results open a way for improvement of 
tunneling transistor, a new type of graphene-based de- 
vice. We believe that our study will stimulate a further 
investigation of these device applications. 



Appendix: Tunneling Hamiltonian 

Below we describe the tunnel-coupled states in 
G/BN/G structure using the 6-column wave function 
[t/'t, 0, V'b] which is written through the spinors V't.f) cor- 
respond to the Gt and Gfc graphene layers. These layers 
are connected through the spinor (j) described the sin- 
gle BN-layer. Within the tight-binding approach, [l^ 
the eigenstate of energy E is determined by the problem 
written through the 6x6 Hamiltonian 



"-GBN 





E hcBN 
hsN — E hsNG 



''BNG 



E 



V'6 



= 0, (A.l) 



where Hq^ and Hq^ are the 2x2 Hamiltonians of Gt, 
Gfc, and BN-layers, while hcBN and Hbng describe weak 
interlayer coupling of Gt and Gf, sheets with BN layer. 
Under a transverse voltage applied, /iCt.t = ± A/2, 
where he is the Hamiltonian of single graphene layer and 
A is a splitting energy between the cross-points in Gt- 
and Gft-layers. For the low-energy {\E\ <1 eV) states, if 
\E\ ^ \£c,v\ where Sc ^SA eV and ~-1.4 eV are the 
c- and w-band extrema energies in BN, we can eliminate 
the spinor (j) from the system (A.l). As a result, the 
eigenstate problem is determined by the 4x4 effective 
tunneling Hamiltonian 



^G + Tt 

T+ 



E 



he 



T 



E 



= 0. 



(A.2) 



Here the 2x2 matrix f describes tunneling through BN 
insulator, while ft and % correspond to the tunneling 
renormalization of t- and 6-states: 



T W —hcBMh^^^hsNG 

Tt ~ —hGBNhg^jhQgjq, h ~ -h 



(A.3) 



BNG^BN^BNG- 



Thus, we arrive to the Hamiltonian (1) with the renor- 
malization contributions Tt.t included to hc^ ^ ; these cor- 
rections are negligible for the weak tunneling regime. 

For the case of A''-layer BN insulator, we consider 
the 2 (A'' -I- 2)-column wave function with A-spinors 
. . . ,(j)N described the BN layers. These states are 
coupled by the interlayer hopping matrices Hbnbn and 
h^NBN which are placed to upper and lower sub- 
diagonals of 2(A + 2) X 2(A + 2) Hamiltonian. ^ Af- 
ter eliminations of the spinors 4>i^ . . . ,(j)N from the tight- 
binding eigenstate problem, we arrive to Eq. (A.3) where 
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the non-diagonal matrix f is replaced by 

TN ~ —hGBNih^]^hBNBN)^~^h^]<fhBNG- (-^-.4) 

For the numerical estimates of T and Tg given by Eq. 
(8) we assume that the diagonal matrix h^]^ is of the 



order of e^]^ where Ebn determined by £c,v The hop- 
ping matrices hcBN, hsNG, and Ubmbn are estimated 
by the interlayer overlap integral 7 and they are strongly 
dependent on a stacking order of G/BN/G structure. 



[1] E. L. Wolf, Principles of Electron Tunneling Spec- 
troscopy, (Oxford University Press Inc., New York, 2012). 

[2] J. Smoliner, E. Gornik, and G. Weiniann, Phys. Rev. 
B 39, 12 937 (1989); J. P. Eisenstein, Superlatt. and 
Microstruct. 12 107 (1992). 

[3] J. Wang, P. H. Bcton, N. Mori, L. Eaves, H. Buhmann, 
L. Mansouri, P. C. Main , T. J. Foster, and M. Henini, 
Phys. Rev. Lett. 73 1146 (1994); N. Mori, P. H. Bcton, 
J. Wang, and L. Eaves, Pliys. Rev. B 51 1735 (1995). 

[4] R. F. Kazarinov and R. A. Suris, Sov. Physics - Semicon- 
ductors, 6, (1972). 

[5] J. P. Eisenstein, T. J. Gramila, L. N. Pfeiffer, and K. W. 
West, Phys. Rev. B 44, 6511 (1991); J. A. Simmons, S. 
K. Lyo, J. F. Klem, M. E. Sherwin, and J. R. Wcndt, 
Phys. Rev. B 47, 15741 (1993). 

[6] L. Zheng and A. H. MacDonald, Phys. Rev. B 47, 10619 
(1993); O. E. Raichev and F. T. Vasko, Journ. of Phys. 
Cond. Matter, 8, 1041 (1996). 

[7] L. A. Ponomarenko, A. K. Geim, A. A. Zhukov, R. Jalil, 
S. V. Morozov, K. S. Novosclov, I. V. Grigoricva, E. 
H. Hill, V.V. Cheianov, V. I. Falko, K. Watanabe, T. 
Taniguchi, and R. V. Gorbachev, Nature Physics 7, 958 
(2011). 

[8] L. Britnell, R. V. Gorbachev, R. Jalil, B. D. Belle, F. 

Schedin, M. I. Katsnelson, L. Eaves, S. V. Morozov, N. 

M. R. Peres, J. Leist, A. K. Geim, K. S. Novoselov, and 

L. A. Ponomarenko, Science 335, 947 (2012). 
[9] K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, 

Y. Zhang, S. V. Dubonos, I. V. Grigorieva, and A. A. 

Firsov, Science 306, 666 (2004). 
[10] S. B. Kumar, G. Seol, and J. Guo, Appl. Phys. Lett. 101, 

033503 (2012). 

[11] For general consideration of the tunneling Hamiltonian 
approach see Ref. 1 and monograph by F. T. Vasko and 



A. Kuznctsov, Electronic States and Optical Transitions 
in Semiconductor Heterostructures, (Springer, New York, 
1998). 

[12] F. T. Vasko and O. E. Raichev, Quantum Kinetic Theory 
and Applications (Springer, New York, 2005). 

[13] T. Ando, J. Phys. Soc. Jpn. 75, 074716 (2006); T. 
Stauber, N. M. R. Peres, and A. H. Castro Neto, Phys. 
Rev. B 78, 085418 (2008); F. T. Vasko and V. V. Mitin, 
Phys. Rev. B 84, 155445 (2011). 

[14] We use the model of disorder with random potentials Vx. 
described by the Gaussian correlator (1414'), with the 
averaged amplitude V, and the correlation length Ic, if 
the cut-off energy Ec — vh/lc exceeds the energy scale 
under consideration. For the low-energy region, \E\ <C 
Ec, the model gives oc E relaxation rate, i. e. we deal 
with a short-range scattering and the coupling constant 
is given hy g — {V Ic/TivYtx /2. 

[15] J. Slawinska, I. Zasada, and Z. Klusek, Phys. Rev. B 81, 
155433 (2010); R. M. Ribeiro and N. M. R. Peres, Phys. 
Rev. B 83, 235312 (2011). 

[16] F. T. Vasko and V. Ryzhii, Phys. Rev. B 76, 233404 
(2007); N. M. R. Peres, Rev. Mod. Phys. 82, 2673 (2010). 

[17] C. R. Dean, A. F. Young, 1. Meric, C. Lee, L. Wang, S. 
Sorgenfrei, K. Watanabe, T. Taniguchi, P. Kim, K. L. 
Shcpard, J. Hone, Nature Nanotechn. 5, 722 (2010). 

[18] From Eq. (9) one obtains /-a = —I a, I-Sep = —Isep, 
and I-ep = lep- According to Eq. (11), A and £Ft,b are 
even functions of voltages, so that I remains the same 
under Vtb,bg -> —Vtb,bg or Vtg,tb,bg ~^ —Vtg,tb,bg- 

[19] F. Amet, J. R. WiUiams, A. G. F. Garcia, M. Yankowitz, 
K. Watanabe, T. Taniguchi, and D. Goldhaber-Gordon, 
Phys. Rev. B 85, 073405 (2012). 



